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THE  EVAPORATION  OF  SESSILE  DROPLETS 
I .  THEORETICAL  FACTORS  FOR  A  TWO-COMPONENT  SYSTEM 


1.  INTRODUCTION 


There  are  many  situations  in  which  evaporating  sessile 
dropleub  contain  more  than  one  component  (e.g.,  agricultural 
spraying,  forestry  spraying,  and  chemical  agent  dissemination). 
As  a  result,  a  model  for  the  evaporation  of  sessile  droplets 
composed  of  two-component  solutions  is  provided  in  this  report. 
The  purpose  of  this  report  is  to  provide  a  model  and  computer 
program  that  can  be  used  for  the  evaporation  of  single-component 
sessile  droplets  in  an  unclassified  and  easy-to-follow  form.  A 
sample  application  to  a  single-component  sessile  droplet  is 
provided.  Also,  the  model  and  computer  program  provided  are 
designed  to  be  used  to  calculate  the  evaporation  of  sessile 
droplets  composed  of  a  two-component  solution.  This  model  was 
developed  for  solving  numerous  mission  problems. 

There  are  many  reports  that  deal  with  related  topics 
(e.g.,  the  evaporation  of  unitary  droplets  into  an  atmosphere 
composed  of  mixtures  of  vapors  and/or  gases),  but  nothing  has 
been  found  that  deals  directly  with  the  topic  discussed  in  this 
report.  Walker  and  Penski^  report  a  model  for  the  evaporation 
of  mixtures  in  a  suspended  droplet.  Coutant  and  Penski^  have 
published  a  report  on  the  evaporation  of  neat  sessile  droplets 
as  a  function  of  wind  velocity  and  contact  angle.  Figure  1 
defines  the  contact  angle  and  illustrates  how  the  Coutant  and 
Penski  model  regard  ic  as  a  constant.  Figure  2  illustrates  a 
decreasing  contact  angle  similar  to  that  considered  in  this 
report.  It  is  intended  that  this  report's  model  will  become  a 
module  in  large  computerized  models.  These  larger  models  will 
consider  additional  situations  including  diffusion  of  the 
droplet  fluid  into  the  surface  as  illustrated  in  Figure  3. 


2  .  FUNDAMENTAL  EQUATIONS 

2.1  Hildebrand  Parameter. 


The  Hildebrand  Parameter,^  5,  is  commonly  called  the 
"solubility  parameter"  and  less  frequently  the  "total  cohesion 
parameter."  It  is  defined  as  follows; 


8  =  [  ( AHv  -  RT)/vl^^^ 


(1) 


where 

AHv  =  molar  enthalpy  of  vaporization 
R  =  ideal  gas  law  constant 
T  =  temperature  (Kelvin) 

V  =  molar  volume  of  liquid 
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Figure  2.  Sessile  Droplet  Evaporation  with  Decreasing  Contact 
Angles:  Hj^  >  B;  >  H3,  t]  <  t2  <  t3. 
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Solut ions . 

According  to  Hildebrand,  Prausnitz,  and  Scott, 3  the 
equations  for  regular  solutions  can  be  written  as  follows: 


2 

RT  In  =  vi  .4>  2  ~  ^2)^  2  A12  ^2] 


(2) 


2  2 

RT  In  =  V2  .4  j  [  (bj  -  {>2^  +  2  Aj2  Sj  82!  (3) 

whe  r  e 

X,  =  activity  coefficient  of  component  i 

i 

Vj^  =  molar  volume  of  component  i 

;  =  volume  fraction  or  volume  of  component  i  divided 

‘  by  total  volume  of  mixture 

6.^  =  Hildebrand  parameter  of  component  i 

A  =  binary  coefficient  of  components  i  and  j 

If  equations  2  and  3  are  solved  for  the  activity 
coefficients,  the  following  equations  result: 


2 

vi  ^.2  f 

Xi  "  ® 


62)"  +  2  A^2 


(4) 


2 

V2  1 j  [ 

X;,  e 


52)2  +  2  A, 2  61  52]/RT 


(5) 


2 . 3  Picknett  and  Bexon  Equations. 

The  Picknett  and  Bexon  equations^  were  developed  for  the 
evaporation  of  a  sessile  droplet  in  a  static  atmosphere.  They 
are  as  fol lows ; 


['  =  0  .  5  K  ETp  (  6  ) 

whe  r  e 

K  -  4  ri  J  (C  -  C„)  (7) 
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D 


diffusion  coefficient  of  vapor  in  atmosphere 


C  =  concentration  of  vapor  in  equilibrium  with  droplet 
fluid  at  surface  of  droplet 

Coo  =  concent  rat  ions  of  vapor  at  a  large  distance  from 

droplet,  usually  very  near  to  zero  in  practical  cases. 

P  =  droplt-t.  dv^nsity 

n  =  3.1416 

CO  =  n  -  ti  -  contact  angle, Brit.  definition  (0) 

6  =  contact  angle  of  droplet  with  surface  (Figure  1) 

E  =  [n/3(l  -  cos  B)‘(2  cos  (9) 

The  term  T  is  defined  as  follows: 

For  0  <  H  <  0.175  radian, 

0.6366  0  t  0.09591  U  '  -  0 . 06 1 4  4  H (10) 

For  0.175  <  B  <  n  radian, 

T=  8.957  10'  +0.6333n+0. 1160  0  •'-0.088780 ’+0.01033  11'  (II) 

2 . 4  Coutant  and  Penski  Equation . 

The  Coutant  and  Penski^  equation  relates  the  evaporation 
of  a  sessile  droplet  to  its  mass,  I',  degree  of  ventilacir;n  (as 
measured  by  its  Reynolds  number),  and  droplet  height.  .•'i  '^u.mber 
of  additional  factors  are  included  in  P. 

-dm/dt  =  I’m*''  [1  +  cRe'’'  (a/h)  '']  (12) 

where 


dm/dt  - 

rate  of 

droplet 

evapor at  ion 

m  = 

d  r  c  p  ’ 

mass 

t  = 

■:  :  me 

c  = 

cons  tan 

t  of  C o u 

tant  and  Penski  equation 

Re  -- 

Reynold 

s  number 

-  2  h  'P/v 

a  =■ 

height  of 

droplet 

h  - 

h'dqh,  t  of 

tunnel 

qj  : 

'Wind  ve 
d  r  '.jp  i  F  t 

1  oc  1  t  y  1 
Lo  tha! 

:  relative  velocity  of  t  ;u 
■  j  f  '.he  at  .mo  s  p  h  e  r  e 

V  ~ 

k 1 nema  t 

1C  'yiscc 

0  i  t  y  .of  the  ..i  t. mo ,s p h e  r  e 

2.5 


Height  of  Sessile  Droplet. 

The  droplet  height^  is  given  by  the  following  equations: 

\  =  1  +  0.75  [2  sin  0/(1  -  cos  0)]^  (14) 

a  =  [  6  m/(  n  p  X)  ( 15  ) 

The  variable  X  is  introduced  only  to  simplify  the  equation  for 
droplet  height. 

2 . 6  Concentration  of  Vapor. 

The  concentration  of  vapor  from  the  evaporating  droplet 
can  be  determined  fairly  accurately  near  ambient  temperatures 
for  the  ideal  gas  law:5 

p°  V  =  nRT  ( 16 } 

where 


=  vapor  pressure  of  fluid 

V  =  volume  of  vapor  at  P° 

n  =  number  of  moles  of  vapor  in  V 

At  the  low  partial  vapor  pressures  under  consideration, 
corrections  for  nonideality  usually  prove  unnecessary.  Equation 
16  can  easily  be  transformed  to  provide  concentration  as 
follows : 


C  =  nM/V  =  MP°/(RT)  (  17 ) 

where  M  =  molecular  weight  of  vapor. 

2 , 7  Partial  Vapor  Pressures. 

The  partial  vapor  pressure  of  component  i  of  a  mixture 
is  given  by  the  following  equation:^ 

o 

P|  =  P  ^  Xj  4>  j  (  i8) 

whe  r  e 

p^  =  partial  pressure  of  component  i 

o 

p  ^  =  vapor  pressure  of  component  i 

=  activity  coefficient  of  component  i 

'  =  volume  fraction  or  volume  of  component  i  divided 

by  total  volume  of  mixture 
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Volume  fractions  are  more  widely  used  for  nonelectro¬ 
lytes  than  mole  fractions  because  volume  fractions  deal  more 
adequately  with  molecular  size  effects. 


3.1.1 


MATHEMATICAL  MODEL 


Starting  Conditions  and  Properties. 
Atmospheric  Conditions  and  Properties, 


3.1.2 


3.1.3 


V  =  kinematic  viscosity  of  the  atmosphere 

=  wind  velocity  or  relative  velocity  of  the  droplet 
to  that  of  the  atmosphere 

T  =  temperature  in  degrees  Kelvin.  In  this  case, 
evaporation  is  assumed  to  be  slow.  Therefore, 
atmospheric  and  droplet  temperature  are  equal. 

C„^  =  concentration  of  vapor  at  a  large  distance  from 
the  droplet  for  each  component,  usually  equal  to 
zero  in  practical  cases. 

h  =  height  of  tunnel 
Droplets . 

mo  =  initial  droplet  mass 
$  ^  =  volume  fraction  of  component  1 
0^  =  volume  fraction  of  component  2 
H  =  contact  angle  of  droplet  with  surface 
Properties  of  Two  Components. 


D  =  diffusion  coefficient  of  component  i  vapor  in 
^  atmosphere 

M^  =  molecular  weight  of  component  i 


Pj  =  vapor  pressure  of  component  i 

8^  =  Hildebrand  parameter  of  component  i 

A^j  =  binary  coefficient  of  components  i  and  j 

=  density  of  component  i 

3  .  2  Computational  System  Overview . 

Calculate  all  the  dependent  variables  by  stepping  from 
Sections  3.2.1  through  3.2.9  in  the  order  given. 
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3.2.1 


Increment  Time  by  a  Small  Number. 


The  calculation  follows  the  procedure  previously 
described  by  Penski^ 

t  -  t.  +  I  (19) 

where 

tj  =  time  after  previous  increment 
T  =  small  increment  in  time 

3.2.2  Activity  Coefficients. 

Calculate  the  activity  coefficients  (equations  4  and  5) 
as  follows: 

v^  $2  ■  ^2)^  +  2  8^  52]/RT 

y^=e  ( 20 ) 

2  7 

V2  (D^  [  (85  -  62)'^  +  2  Ai2  Sj  821/RT 

72^^  ( 21 ) 

3.2.3  Drop  Volume. 

Q  =  "'/Pi 

where 

Q  =  volume  of  droplet 

=  density  of  component  i 

3.2.4  Reynolds  Number. 

Re  =  2  h  ^?/v  (23) 

3.2.5  Concentration  of  Vapor  at  Droplet  Surface . 

The  partial  vapor  pressure  of  component  i  of  a  mixture 
is  given  by  the  following  equation: 


’  y  i) 

1  '  \ 


and  the  concentration  of  vapor  of  the  two  components  is  as 
follows : 


(24) 

(25) 
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3.2.6 


where 


Cj  =  Ml  Pi/(RT)  (26) 

=  M2  P2/tRT)  (27) 

Cl  =  concentration  of  vapor  of  component  i  in 

equilibrium  with  the  fluid  at  surface  of  droplet. 

M|  =  molecular  weight  of  vapor  of  component  i 


Fj  =  0.5  KjE  Tp 
r2  =  0.5  K2E  Tp 


(28) 

(29) 


Kj  =  4  nDi(Ci  -  C„i)  (30) 

K2  =  4  nD2(C2  -  C„2)  (31) 

CO  =  n  -  H  (32) 

E  =  [n/3(l  -  cos  H)2  (2  +  cos  B)]-^/^  (33) 

For  0  <  B  <  0.175  radian, 

T  =  0.6366  B  +  0.09591  B^  -  0.06144  B^  (34) 

For  0.175  <  6  <  t  radian, 


T=  8.957  10' ^  +  0.63  33;'.  +  0.11 60  0^-0. 08878  (35) 

3.2.7  Height  of  Sessile  Droplet. 

The  droplet  height2  is  given  by  the  following  equations: 
X  =  I  +  0.75  [2  sin  B/(l  -  cos  B)]^  (36) 

a  =  [  6  m/(n  p  X)  (  37  ) 

3.2.8  Mass  of  Droplet  and  Mass  of  Each  Component. 

The  Coutant  and  Penski  equation  relates  the  evaporation 
of  a  sessile  droplet  to  its  mass,  Tj,  and  degree  of  ventilation 
as  measured  by  its  Reynolds  number: 

mi  =  mij  -  Ti  m^/^(l  +  c  Re^'-'^  (a/h)'^‘‘’j  t  (38) 
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"^2  "  "'2j  ~  ^2  ^  (a/h)°*^]  x 

(39) 

m  =  m^  +  m2 

where  m 

ij  =  mass  of  component  i  for  last  iteration. 

3.2.9 

Volume  Fraction. 

$1  ==  nij/(Pj  Q) 

(41) 

^2  =  m2/(p2 

(42) 

3.2.10 

Contact  Angle  Decay. 

Provision  for  contact  angle  decay  is  provided  in  the 
model.  The  following  equation  is  used. 


B  Bq(Ao  -  Aj  m/mo)  (43) 

where  Bq  is  the  initial  contact  angle,  and  Aq  and  Ai  are 
constants . 

At  this  point,  calculation  returns  to  Section  3.2.1  and 
starts  over  providing  all  the  variables  as  a  function  of  time 
through  multiple  repetitions. 


4 .  RESULTS 

The  definitions  of  the  terms  used  in  the  program  and 
their  units  are  provided  in  Appendix  A.  Appendix  B  provides  the 
computer  program  developed  in  this  effort.  Figure  4  compares 
data  collected  by  experimental  methodology  by  the  authors®  with 
calculations  performed  with  the  model  and  computer  program 
described  in  this  report.  The  data  used  was  for  a  sessile  water 
droplet  on  an  aluminum  surface  in  a  small  wind  tunnel  with  a  low 
air  flow.  Details  of  the  calculations  are  provided  in  Appendix 
C. 


5.  DISCUSSION 

To  keep  this  report  unclassified  and  relatively  simple, 
only  water  data  is  included.  By  taking  the  concentrations  of 
the  second  component  as  zero,  the  model  may  be  applied  to  neat 
liquids . 


If  equations  are  introduced  to  vary  the  properties  of 
the  components,  the  model  may  be  used  for  conditions  other  than 

isothermal.  For  example,  the  Antoine  equation^  may  be  used  to 
o 

vary  p^  and  (vapor  pressures  and  enthalpy  of  vaporization 

for  each  component). 
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6. 


CONCLUSIONS 


This  report  provides  a  model  and  computer  program  that 
can  be  used  to  calculate  evaporation  of  single-component  sessile 
droplets . 


A  model  and  computer  program  are  provided  that  can  be 
used  to  calculate  the  evaporation  of  sessile  droplets  composed 
of  two-component  solutions. 

Data  for  the  evaporation  of  a  sessile  water  droplet  on 
an  aluminum  surface  in  a  small  wind  tunnel  with  a  low  air  flow 
compares  favorably  with  calculations  using  the  model  and 
computer  program. 
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GLOSSARY 


a  =  Height  of  droplet 

C-.  =  Concentration  of  vapor  of  component  i  in  equilibrium 

with  the  fluid  at  surface  of  droplet 

=  Concentration  of  vapor  at  a  large  distance  from 

droplet  for  each  component,  usually  equal  to  zero  in 
practical  cases 

c  =  constant  of  Coutant  and  Penski  equation 
=  molar  enthalpy  of  vaporization 
D  =  diffusion  coefficient  of  vapor  in  atmosphere 
E  =  defined  by  equation  9 
h  -  height  of  tunnel 
K  -  defined  by  equation  7 
M;  =  molecular  weight  of  component  i 
m  -  droplet  mass 

m,.  -  initial  droplet  mass 

m|  =  mass  of  component  i  in  droplet 

m. ,  j  =  mass  of  component  i  for  last  iteration 

n  =  number  of  moles  of  vapor 


O 

Pj  =  vapor  pressure  of  component  i 
Pj  =  partial  pressure  of  component  i 
R  =  ideal  gas  law  constant 
Re  =  Reynolds  number  -  2  h  'V/v 
T  =  teniperature  in  Kelvin 
t  =  time 

t.  =  time  after  last  increment 

j 
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=  molar  volume  of  component  i 

V  =  volume  of  vapor 
GREEK  TERMS 

8^  =  Hildebrand  parameter  of  component  i 

<fi  =  volume  fraction  or  volume  of  component  i  divided  by 
total  volume  of  mixture 

r  =  defined  by  equation  ft 

y  =  activity  coefficient  cf  component  i 

0  =  contact  angle  of  droplet  with  surface 

9o  =  initial  contact  angle  of  droplet  with  surface 

A =  binary  coefficient  of  components  i  and  ] 

X  =  defined  by  equation  14 

V  =  kinematic  viscosity  of  the  atmosphere 
n  =  3.1416 

'F  =  wind  velocity  or  relative  velocity  of  the  droplet  t 
that  of  the  atmosphere 

=  density  of  component  i 

T  =  small  increment  in  time 

Q  =  volume  of  droplet 

T=  defined  by  equations  10  and  11 


24 


APPENDIX  A 

DEFINITIONS  OF  COMPUTER  TERMS 


25 


Blank 


26 


A1  Height  of  droplet,  cm 

Cl  Concentration  of  component  1  in  vapor,  g/cm^ 

C2  Concentration  of  component  2  in  vapor,  g/cm^ 

C4  Concentration  of  component  1  in  atmosphere  leaving 

tunnel,  ug/cm  ^ 

C5  Concentration  of  component  2  in  atmosphere  leaving 

tunnel,  yg/cm  ^ 

C6  Constant  of  Coutant  and  Penski  equation 

D1  Density  of  component  1,  g/cm^ 

D2  Density  of  component  2,  g/cm^ 

D5  Density  of  atmosphere,  g/cm^ 

D8  Diffusion  coefficient  of  component  1  in  atmosphere, 

cm^s 

D9  Diffusion  coefficient  of  component  2  in  atmosphere, 

cm2  3  - 1 

E8  T,  see  equations  10  and  11 

E9  E,  see  equation  9 

FI  ,  see  equation  28 

F2  see  equation  29 

F3  Right  most  two  factors  in  equation  12 

F4  r^/Ki 

G3  Activity  coefficient  of  component  1,  unitless 

G4  Activity  coefficient  of  component  2,  unitless 

HI  Height  of  tunnel,  cm 

12  Number  of  points  per  data  set 

13  Point  counter 

17  Data  set  number 

K  273.16 
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K1  Temperature  of  droplet,  Kelvin 

K2  Atmospheric  temperature,  Kelvin 

K5  Kj,  equation  30,  g/(cm  s) 

K6  K2 ,  equation  31,  g/(cm  s) 

L  Binary  coefficient,  unitless 

M  Droplet  mass,  experimental,  mg  as  input  converted  to  g 

Ml  Mass  of  component  1,  calculated,  g 

M2  Mass  of  component  2,  calculated,  g 

M4  Molecular  weight  of  component  1,  g 

M5  Initial  droplet  mass,  g 

M7  Droplet  mass  for  previous  iteration,  g 

M8  Molecular  weight  cl  component  2,  g 

N  Viscosity  of  atmosphere,  poise 

O  Initial  volume  of  droplet  minus  volume  of  impurities, 

cm^ 

04  Initial  volume  of  component  1,  cm^ 

05  Initial  volume  of  component  5,  cm^ 

P  Total  of  component  partial  vapor  pressures,  Torr 

PI  Partial  vapor  pressure  of  component  1,  Torr 

P2  Partial  vapor  pressure  of  component  2,  Torr 

P3  Pressure  of  atmosphere,  Torr 

P4  Purity  of  component  1  in  weight  fraction,  unit  less 

P5  Purity  of  component  2  in  weight  fraction,  unitless 

P7  Vapor  pressure  of  component  1,  Torr 

P9  Vapor  pressure  of  component  2,  Torr 

R  Ideal  gas  law  constant,  cm^  Torr/(Kelvin  mole^ 

R2  Ideal  gas  law  constant,  cal/(KeIvin  mole) 
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R3 

Reynolds  number,  unitless 

S8 

Hildebrand  parameter  of  component  1, 

(cal/cm^) 

1/2 

S9 

Hildebrand  parameter  of  component  2, 

( cal/cm-^) 

i/2 

T 

Measured  time,  input  as  min,  converted  and  used  as  s 

T1 

Droplet  temperature,  C 

T2 

Atmospheric  temperature,  C 

T6 

Coiii-dct  angle,  inli.ially  iii  deg  then 

ccnvcr ted 

1  ^  r  aH 

T7 

Initial  contact  angle,  deg 

T8 

Calculated  time,  s 

T9 

Increment  in  time,  s 

V 

Kinematic  viscosity  of  atmosphere,  cm^/s 

VI 

Atmospheric  velocity,  cm/s 

V9 

Droplet  volume,  cm^ 

V4 

Volume  of  component  1,  cm^ 

V5 

Volume  of  component  2,  cm^ 

W1 

Increment  of  mass  of  component  1,  g 

W2 

Increment  of  mass  of  component  2,  g 

W3 

Droplet  mass,  g 

XI 

Volume  fraction  of  component  1,  g 

X2 

Volume  fraction  of  component  2,  g 

X4 

Weight  of  component  1  divided  by  initial  mass 
droplet,  unitless 

of 

X5 

Weight  of  component  2  divided  by  initial  mass 
droplet,  unitless 

of 

Z 

Concentration  of  vapor  component  in 
tunnel  per  gram  of  component  leaving 

atmosphe  r  c 
droplet , 

■  leaving 
mg/cm  V9 
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2  .  2  Activity  Coefficients  of  Regular  Two-Component 

Solutions ■ 

According  to  Hildebrand,  Prausnitz,  and  Scott, 3  the 
equations  for  regular  solutions  can  be  written  as  follows: 


RT  In  =  vi  $  2  [(8^  -  62)2  +  2  A12  83] 


RT  in  Y  =  V2  [(8j  -  82)^  +  2  A^2  83] 


where 


/.  =  activity  coefficient  of  component  i 

=  molar  volume  of  component  i 

$  .  =  volume  fraction  or  volume  of  component  i  divided 
^  by  total  volume  of  mixture 

8j^  =  Hildebrand  parameter  of  component  i 

Aj^j  =  binary  coefficient  of  components  i  and  j 

If  equations  2  and  3  are  solved  for  the  activity 
coefficients,  the  following  equations  result: 


VI  $2  n8i  -  82)  +  2  A^2  Si  52]/RT 


r,  =  e 


V2  $ ,  [  (8^  -  82)2  +  2  Ai2  81  82]/RT 

~  /a  ^ 


2 . 3  Picknett  and  Bexon  Equations. 

The  Picknett  and  Bexon  equations^  were  developed  for  the 
evaporation  of  a  sessile  droplet  in  a  static  atmosphere.  They 
are  as  follows: 


r  =  0,5  K  ETp 


where 


K  =  4  n  D  (C  -  CJ 


11 


D  =  diffusion  coefficient  of  vapor  in  atmosphere 

C  =  concentration  of  vapor  in  equilibrium  with  droplet 
fluid  at  surface  of  droplet 

C«,  =  concentrations  of  vapor  at  a  large  distance  from 

droplet,  usually  very  near  to  zero  in  practical  cases. 


P  =  droplet  density 
n  =  3.1416 

CO  =  n-9  =  contact  angle, British  definition  (8) 

6  =  contact  angle  e®  droplet  with  surface  (Figure  1) 

E  =  [n/3(l  -  cos  0)^(2  +  cos  6)]"^'^^  (9) 

The  termT  is  defined  as  follows: 

For  0  <  0  <  0.175  radian, 

T=  0.6366  6  +  0.0S591  0^  -  0.06144  0^  (10) 

For  0.175  <  0  <  n  radian, 

T=  8.957  10‘5  +  0.6333  e  +  0.1160  0^  -  0.08878  +  0.01033  o'*  (11) 

2 . 4  Coutant  and  Penski  Equation. 


The  Coutant  and  Penski2  equation  relates  the  evaporation 
of  a  sessile  droplet  to  its  mass,  F,  degree  of  ventilation  (as 
measured  by  its  Reynolds  number),  and  droplet  height.  A  number 
of  additional  factors  are  included  in  F. 

-dm/dt  =  Fm^^^Il  +  cRe®'®^  (a/h)°-^]  (12) 

where 

-dm/dt  =  ■  ite  of  droplet  evaporation 
m  =  oroplet  mass 
t  =  time 

c  =  constant  of  Coutant  and  Penski  equation 
Re  =  Reynolds  number  =  2  h  'F/u  (13) 

a  =  height  of  droplet 
h  =  height  of  tunnel 

'F  =  wind  velocity  or  relative  velocity  of  the 
droplet  to  that  of  the  atmosphere 

V  =  kinematic  viscosity  of  the  atmosphere 


12 


3.2.1 

Increnienc  Time  by  a  Small  Number. 

The  calculation  follows  the  procedure  previously 
described  by  Penski^ 

t  =  tj  .  r 

(19) 

where 

tj  =  time  after  previous  increment 

T  =  small  increment  in  time 

3.2.2 

Activity  Coefficients. 

Calculate  the  activity  coefficients  (equations  4 
as  follows; 

and  5 ) 

Vi  $2  t  2  Ai2  8^  S^J/RT 

Xi  =  e 

(20) 

V2  $1  [  (8i  -  82)^  +  2  Ai2  8j  S^l/RT 

72  =  e 

(21) 

3.2.3 

Drop  Volume. 

a  =  mj/pj  +  m^/p^ 

(22) 

where 

Q  =  volume  of  droplet 

p^  =  density  of  component  i 

3.2,4 

Reynolds  Number. 

Re  =  2  h  W/v 

(23) 

3.2.5 

Concentration  of  Vapor  at  Droplet  Surface. 

The  partial  vapor  pressure  of  component  i  of  a  mixture 
is  given  by  the  following  equation; 


Pj  “  Pi  ’'i  «i 


and  the  concentration  of  vapor  of  the  two  components  is  as 
follows : 


(24) 

(25) 
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-  Mj  P;/(RT)  (26) 

C2  =  M2  P2/{RT)  (27) 

=  concentration  of  vapor  of  component  i  in 

equilibrium  with  the  fluid  at  surface  of  droplet. 

Mj^  =  molecular  weight  of  vapor  of  component  i 

3.2.6  Picknett  and  Bexon  Equations. 

=  0.5  K^E  Tp""'"'’  (28) 

r,  =  0.5  K.E  (29) 


where 


=  4  nDi(Ci  -  C»i)  (30) 

K2  =  4  nD2(C2  -  C„2)  (31) 

CO  =  n  -  0  (  32  ) 

E  [n/3(l  -  cos  0)’  (2  +  cos  0)]-^/^  (33) 

For  0  <  0  <  0.175  radian, 

T  =  0.6366  0  +•  0.09591  0^  -  0.06144  0^  (34) 

For  0.175  <  0  <  t  radian, 

T  =  8.957  10"^  +  0.6333  8  +  0.1160  0^  -  0.08878  8^+  0.01033  8'*  (35) 

3.2.7  Height  of  Sessile  Droplet. 

The  droplet  height2  is  given  by  the  following  equations; 

\  =  1  +  0.75  [2  sin  0/(1  -  cos  0)]^  (36) 

a  =  [6  m/(npA)]^/3  (37) 

3.2.8  Mass  of  Droplet  and  Mass  of  Each  Component. 

The  Coutant  and  Penski  equation  relates  the  evaporation 
of  a  sessile  droplet  to  its  mass,  Fi,  and  degree  of  ventilation 
as  measured  by  its  Reynolds  number; 

~  ~  c  Re°‘^^  (a/h)°'^]  X  (38) 
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